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Abstract — We present the derivation of post-processing SNR 
for Minimum-Mean-Squared-Error (MMSE) receivers witli im- 
perfect channel estimates, and show that it is an accurate 
indicator of the error rate performance of MIMO systems in the 
presence of channel estimation error. Simulation results show the 
tightness of the analysis. 

Index Terms — MIMO, MMSE receiver, post-processing SNR. 



I. Introduction 

THE key component of a Multiple-Input-Multiple-Output 
(MIMO) communication system in terms of performance 
and complexity is the MIMO detector, which is used for 
separating independent data streams at the receiver. The Maxi- 
mum Likelihood (ML) detectors achieve the optimal error rate 
performance. However, these types of detectors, including the 
near-optimal sphere decoder and its variants, are usually not 
suitable for practical systems due to their high complexity. Lin- 
ear detectors, such as Zero-Forcing (ZF) and MMSE, achieve 
suboptimal performance, however, they are widely used in 
practical systems due to their low complexity implementations. 
Among linear receivers, MMSE is the optimal solution and 
seems to be the mainstream implementation choice due to its 
superior performance over ZF detectors. 

Perfect channel state information (CSI) is usually assumed 
in the literature when simulating or analyzing the performance 
of linear detectors |[T|, Q. However, in practice the channel 
estimates are inherently noisy. Important work Q, ||4l has 
characterized the error rate performance of ZF receivers in 
the presence of channel estimation error Nevertheless, less is 
known for the case of MMSE detectors in practical scenarios. 
For ZF and MMSE receivers, the joint effect of phase noise 
and channel estimation error is considered in |5| and the 
performance is analyzed in terms of the degradation in signal- 
to-noise-plus-interference-ratio (SINR) without expressing the 
closed form performance indicators or error rate analysis. The 
SINR derivations for the MMSE case in |5| are done only for 
low SNR region. In both |3| and |5|, channel estimation error 
variance is assumed to be constant for all SNRs. This is not 
realistic approach for packet based or bursty communication 
systems as the channel estimation error is in fact a function of 
the SNR. In this letter, we analyze the MMSE receivers in the 
presence of channel estimation error, and derive a closed form 
post-processing SNR expression, which provides an accurate 
estimate of the error rate performance. The error rate perfor- 
mance is investigated for both the constant channel estimation 
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error variance case and the case with a realistic channel 
estimation algorithm where the estimation error variance is 
clearly dependent on the channel SNR. We believe that it is 
a very useful tool for throughput prediction in link adaptation 
protocols and for error rate analysis in general. Accuracy of 
the analytical results is verified through simulations. 

II. System Model and Derivations 

We consider a MIMO system where the transmitter is 
equipped with Nt antennas, and the receiver uses Nr antennas. 
The Nr y. \ received signal vector y can be expressed as 



Hx 



(1) 



where x is the transmitted signal vector, H is the x Nt 
channel matrix, and n is the Nj. y- 1 additive Gaussian noise 
vector with zero mean and covariance matrix E [nn^] = 
NqI. We assume an uncorrelated Rayleigh flat channel, i.e. 
entries of H are i.i.d. zero mean circularly symmetric complex 
Gaussians (ZMCSCG) with unit variance, and the signal 
energy at each transmit antenna is assumed to be equal to 
Es- 

The receiver can estimate the transmitted signal vector by 
applying the MMSE detector to the received signal, x = 
Wy = WHx + Wn. Using the orthogonality principle [6J, 
the MMSE detector W is derived as 
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At the output of the MMSE detector, the residual signal plus 
interference from other spatial streams is well approximated 
as Gaussian [7 1 and the post-processing SNR (PPSNR) of /c*'' 
spatial stream is calculated asM 
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The PPSNR is a good indicator for the error rate perfor- 
mance of MIMO systems, and therefore employed in link 
adaptation algorithms to predict the uncoded error rate l?). 
Since the output of the MMSE detector is Gaussian, the 
bit error rate of a specific modulation can be calculated by 
simply plugging the PPSNR value into the AWGN error rate 
formula of the modulation. The same technique is also used 
for theoretical derivation of error rate performance in fading 
channels. 

' (. . ; denotes the (fc, l)*^ entry of the matrix. 
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This definition of PPSNR holds if the channel is perfectly 
known at the receiver. However, in practice, the channel 
matrix has to be estimated by the receiver, and the estimated 
channel is inherently noisy in practical systems. We model the 
estimated channel matrix as 



H = H + AH 



(4) 



where AH denotes the estimation error matrix which is 
uncorrected with H, and its entries are ZMCSCG with 
variance Cg. The quality of channel estimation is captured 
by (Tg, which can be appropriately estimated depending on 
the channel estimation method. We assume that each block 
(packet), that undergoes a specific channel realization, H, 
observes a different realization of AH at the receiver. This 
situation occurs in packet based communication systems like 
802. lln where the channel is estimated on a per packet basis. 

A. PPSNR derivation for practical systems 

In this section, we derive the PPSNR for practical MIMO 
systems which observe channel estimation error. The receiver 
uses the estimated channel H to calculate the MMSE detector 
as 
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We write the imperfect MMSE solution as W = W + AW. 
Now, the MMSE estimate of the signal vector becomes 

X = (W + AW) y = WHx + A WHx + Wn + A Wn 

signal post-detection noise 

(6) 

We observe that there are additional interference and noise 
terms caused by AW, and denote the post detection noise as 
n = AWHx + Wn + AWn. With this definition for the post 
detection noise, the PPSNR of the fc*'' spatial stream in the 
presence of channel estimation error can be expressed as 
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where we replaced the original noise covariance in ^ with 
the covariance of n, which is calculated as 



E [r 



E [AWHxx^H^AW^] + E [Wnn-^W^] 

+ E [Wnn^AW^] + E [AWnn^W^] 

+ E [AWnn^AW^] (8) 



In order to calculate the terms in (|8]l, we need to first derive 
AW. For small al, the AH^AH term in (|5]l becomes 
negligible compared to others. Hence, we can rewrite (|5| as 
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which can be further simplified using the matrix approximation 
(P + e2Q)~^ p-i-e2p-iQp-i for small e^. Let us also 



define K {yL"U. + fai) for brevity and simplify Q as 

W = [K - K (H-f^AH + AH^H) K] (H + AH)^(10) 
= KH^ K (H^AH + AH^H) KH^ 

+ KAH^ - K (H^AH + AH^H) KAH^ (11) 



small compared to the other terms 

Finally the desired error matrix becomes 

AW = -K(H^AH + AH^H)KH^ + KAH^. (12) 

Using the above approximation, we can now calculate the 
terms in (|8]l. We first note that the third and fourth terms 
in (|8]l are zero since E [AW] ^ 0. The second terms is 
E [Wnn^W-^] = A^qWW-^^, and the first term becomes 
E [AWHxx^H^AW^] = E^E [AWHH^AW^] . Be- 
low, we calculate the first and last terms in ([8]) by plugging 
the error matrix ( [T2] i into (|8]l. 

E [AWHH^AW-^] 

= E [KH^AHKH^HH^HK^AH^HK^] 
+ E [KH^AHKH^HH^HK^H^AHK^] 

- E [KH^AHKH^HH^AHK^] 

+ E [KAH^HKH^HH^HK^AH^HK^] 
+ E [KAH^HKH^HH^HK^H^AHK^] 

- E [KAH^HKH^HH^AHK^] 

- E [KAH^HH^HK^AH^HK^] 

- E [KAH^HH^HK^H^AHK^] 

+ E [KAH^HH^AHK^] (13) 

It can be proven that E [AH A AH] = E [AH^AAH^] = 
for any deterministic matrix A. Hence the second, third, fourth 
and seventh terms in (13) are zero. For the remaining terms 
we use the fact that E [AHAAH^] = a^tr (A) I, and obtain 

E [AWHH^AW^] 



cr,^tr (KH^HH" HK" ) KH" HK 



a'itr ( HKH^HH^ HK" H" ) KK 



-crXHKH"HH")KK^ (14) 
- aftr (HH^HK^H^) KK^ + a^tr (HH^) KK^ 

Similarly, the last term in ([8]), E [AWnn^AW^], can be 
computed following the same way. 

E [AWnn^AW^] = NqE [AWAW^] (15) 
E [AWAW^] 

= cr> (KH^HK-^) KH^HK^ 
+ CT^tr (HKH^HK-f^H^) KK^ 

- CT^tr (HKH^) KK^ 

- CT^tr (HK^H^) KK^ + a^N^-KK" (16) 

Finally, we plug E [nn^] into ^ and obtain the PPSNR in 
the presence of channel estimation error as ( [T7] i. 

The BER of the system in the presence of channel es- 
timation error can be found simply by plugging jk as the 
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1 S.ctXKH^HH^HK^) KH^HK^ + £;^CT2ti-(HKH^HH^HK^H^) KK^ \ 
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(17) 



symbol SNR into the AWGN BER formulas. For example, 
the BER of fc*'' stream for BPSK is = 



and 



Q 



coded 16QAM. 
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III. Results 



In order to test the performance of the analysis, we sim- 
ulated transmission of thousands of packets through uncorre- 
cted Rayleigh flat fading channels. For each SNR point on 
the BER plots, we randomly generate 1000 i.i.d. realizations 
of the channel matrix H. For each specific realization of the 
channel, we transmit 500 packets each of which carries 2000 
information symbols. We perform channel estimation for each 
packet as explained below in Case 1. 

Case 1: In our simulations, we employed the maximum 
likelihood (ML) channel estimation (CE) algorithm, in which 
the channel estimate is obtained via training symbols that are 
known to the receiver. During the training phase, the Nt x Ntr 
training matrix X-tr is transmitted where Ntr > Nt is the 
number of training symbols. The Nr x Ntr received signal is 
Ytr — HXf,, + W where W is the Nr x Ntr noise matrix. 
Then, the ML estimate of the channel is given as |[8| 



H — Ytr^u- i^tr^tr) 



(18) 



It was shown that the optimal training signal has the property 
of XfrX^ = EsNtrl- When this orthogonal training signal is 
employed, the entries of AH are i.i.d. with CN{0,(t'^), and 
the channel estimation noise varianca^is crj 



The estimation error in this case is caused by the AWGN in 
this case. 

The following training signal, which is taken from 802.1 In 
standard [9J, was employed in the simulations. X^^ = \/E^P 
where P is the submatrix formed by first Nt rows and first Ntr 
columns of the bigger matri)|^P, i.e. P = P [1 : Nt;l : Ntr]- 
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It should be noted that with this choice of the training matrix, 
the ML channel estimation at the receiver becomes a very 
simple operation since the matrix inversion, (Xfj-X^) , is 
now a trivial operation. 
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(T„ can also be defined as ai 



Nt 



Nt E /No depending on SNR definition. 
'Ttie P matrix liere is for maximum of 4 spatial streams since tlie standard 
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supports up to 4 streams. Ntr = 4 for Nt > 3. 



Fig. 1. BER for BPSK. Et/No is per tx antenna. 



We present the simulation results for BPSK in Fig. [T] and 
16QAM in Fig. |2]with 1x4, 2x4, 4x4 MIMO configurations. 
Ntr = 4 is used in all the simulations. The case of <t^ — 0, 
i.e. perfect channel estimation, is also included in the results. 

For each channel instance, analytical BER results are ob- 
tained by using the PPSNR derived in the previous section. 
Then, these BERs are averaged over all realizations of the 
channel. 

First thing to notice in Fig. [l] and Fig. |2]is that for = 0, 
simulation and analysis curves exactly match. Performance is 
significantly degraded for the systems experiencing channel 
estimation errors. This is particularly evident for the 4x4 
configurations. 

As it can be seen in Fig. [Tjand Fig. |2j our analysis gives a 
very tight approximation of the real performance. For BPSK 
4x4, and all of the 16QAM configurations the analysis 
results exactly match the simulated performances. For BPSK 
1x4, and 2x4 configurations the analysis results are upper- 
bounds to the real performance at high SNR, however, they are 
still very close to the real performances. The analysis results 
become tighter for higher order modulations and higher order 
MIMO configurations. This is because of the fact that the 
Gaussian assumption, which is made for the post-detection 
noise, is more valid at higher order modulations and MIMO 
configurations. At low SNRs, the total post detection noise n is 
dominated by the additive white Gaussian noise component n 
therefore the assumption is valid even for lower configurations. 
However, at high SNRs the residual interference components 
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Simulation, = 1), perfect CE 

' ' 'O ' ' Analysis, ct^ = 0. perfect CE 1x4-' 
— — — Simulation, ML CE 
Analysis. ML CE 



Ei,/No (clB) 



Fig. 2. BER for 16QAM. Eb/No is per tx antenna. 




0.2 MMSE Simulation 
n.2 MMSE Analysis 
0.2 ZP Simulation 
0.1 MMSE Siimilation 
0.1 MMSE Analysis 
0.1 ZP Simulation 
0.05 MMSE Sinmlntior 
0.05 MMSE Analysis 
0.05 ZP Simulation 
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MMSE Analysis 
ZF Simulation 
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Fig. 3. BER for 4x5 QPSK for different (Te values, = corresponds to 
perfect CE. ZF curves are included for comparison purposes with |3 |. Eb/No 
is per tx antenna. These parameter values are taken from Fig. 2 of t3J for 
comparison purposes. 



from other spatial streams becomes dominant and fi is loosely 
approximated as Gaussian for lower order constellations and 
MIMO configurations. 

It is interesting to note that in contrast to the results obtained 
for ZF detector by [3|, we do not observe any error floor 
on the performance. This is due to the fact that the channel 
estimation error variance for ML estimation gets smaller 
as SNR increases. This is the situation that occurs in practical 
packet based or bursty communication systems where the 
channel estimation is performed for every packet prior to data 
detection, and hence experiences the same noise variance as 
the data transmission. Therefore the channel estimation quality 
is dependent on the SNR. On the other hand, error floors are 
observed in [3| because of the assumption that remains 
constant independent of the SNR. This case is investigated 
below in Case 2. 



Case 2: In addition to ML channel estimation results, we 
also performed simulations with constant erg. Unlike the first 
case, the channel estimation quality is independent of the SNR. 
This situation might arise either when there is a ready channel 
estimate to be used by the receiver formed elsewhere with a 
different additive noise variance, or the channel estimation is 
outdated and the major error in the channel estimation comes 
from the mobility changes in the channel. 

In Fig. [3] the BER performance of a QPSK 4x5 system 
is investigated for CTg — 5%, 10%, 20% using the estimation 
error model in (4). Each packet observes a different realization 
of the random matrix AH with the designated variance a^. 
As expected, we observe error floor in the performance due 
to the constant estimation error variance as in the ZF detector 
case studied in |3|. More importantly, these error floors are the 
same as the ones observed by ZF detector because of the fact 
that the MMSE and ZF detectors exhibit the same behaviour at 
asymptotically high SNR. The simulation results in this case 
also agree with the analysis. 

IV. Conclusion 

In this letter, we presented the analysis of post-processing 
SNR for practical MIMO MMSE receivers which experience 
imperfect channel estimation. Performance of MMSE receivers 
in the presence of channel estimation error is investigated and 
shown to be accurately estimated via analytical results. We 
verified the tightness of the analytical results via simulations. 

Besides the theoretical contributions, we believe that our 
closed form PPSNR expression can be useful for link adap- 
tation purposes in real MIMO systems. There exist link 
adaptation algorithms Q, ifTOl based on PPSNR, however 
perfect CSI is always assumed which might lead to incorrect 
prediction of the throughput. More accurate prediction can be 
achieved using the results presented in this paper 
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